Measuring correlated electron dynamics with time-resolved photoemission 

spectroscopy 
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Time-resolved photoemission experiments can reveal fascinating quantum dynamics of correlated 
electrons. However, the thermalization of the electronic system is typically so fast that very short 
probe pulses are necessary to resolve the time evolution of the quantum state, and this leads to poor 
energy resolution due to the energy-time uncertainty relation. Although the photoemission inten- 
sity can be calculated from the nonequilibrium electronic Green functions, the converse procedure 
is therefore difficult. We analyze a hypothetical time-resolved photoemission experiment on a cor- 
related electronic system, described by the Falicov-Kimball model in dynamical mean-field theory, 
which relaxes between metallic and insulating phases. We find that the real-time Green function 
which describes the transient behavior during the buildup of the metallic state cannot be determined 
directly from the photoemission signal. On the other hand, the characteristic collapse-and-revival 
oscillations of an excited Mott insulator can be observed as oscillating weight in the center of the 
Mott gap in the time-dependent photoemission spectrum. 

PACS numbers: 71.27.+a, 78.47.-p 



I. INTRODUCTION 

Pump-probe experiments with femtosecond time res- 
olution can record various nonequilibrium processes in 
solids directly in the time domain, including those in- 
duced by the Coulomb interaction between electrons, or 
the scattering of electrons on defects and phonons. In 
these experiments excitation of the sample and charac- 
terization of the excited state are accomplished by two 
distinct laser pulses (pump and probe) which hit the 
sample with controlled time delay, and either optical or 
photoemission spectroscopy may be used as probe tech- 
nique. The pump-probe setup has been used to investi- 
gate the dynamics of molecules ji semiconductors^ and 
metals^ for more than two decades. More recently, such 
time-resolved experiments were also performed on several 
strongly correlated materials close to a phase transition, 
where many degrees of freedom contribute to the dynam- 
ics on very different time scales 4^S»L2iS 

In such nonequilibrium solid-state experiments it is a 
major challenge to distinguish the electronic dynamics 
from other degrees of freedom. This is crucial in par- 
ticular for the Mott metal-insulator transition^ which 
is driven by the Coulomb interaction between electrons 
moving in a crystal lattice. An entirely new perspective 
on this phenomenon would open up if one could observe 
the transition as it happens in real time and, e.g., mon- 
itor the formation of well-defined quasiparticles as the 
system goes from an insulating to a metallic state. In 
fact, in several Mott and charge-transfer insulators the 
transition to a metallic state can be induced by a laser 
pump pulse^i 7 -^ So far these experiments have focused 
on the relaxation back to the insulating state, which in- 
volves coupling to degrees of freedom other than the va- 
lence band electrons, and is much slower than the buildup 
of the metallic state after the pump pulse. In the sim- 



plest case the observed relaxation is described by the two- 
temperature models i.e., as the cooling of a hot electron 
gas which is coupled to the colder latticeiiS The true dy- 
namics of the electronic system has so far been observed 
only in simple metals, by looking at the thermalization 
of pump-excited electron distributions due to electron- 
electron scattering 4>H In strongly correlated materials, 
thermalization is apparently much faster. Sufficient time 
resolution is now becoming available due to recent ad- 
vances in femtosecond laser techniques ; 14 ! 15 which have 
already allowed to investigate some solid-state systems 
even on the attosecond time scaled 

In equilibrium, electronic properties of correlated ma- 
terials can be obtained directly from conventional photo- 
emission spectroscopy with continuous light beamsi 17 ! 18 
By contrast, time-resolved measurements are likely to be 
restricted by the frequency-time uncertainty of the probe 
pulse 1 : The energy c of occupied states in the solid from 
which photoclcctrons arc released is determined from the 
kinetic energy of the photoelectrons, the work function 
of the solid, and the photon energy E 7 = ftu] when the 
measurement pulse has finite duration St, the latter is 
determined only up to an uncertainty 8Ej > h/St. In 
a strongly correlated electron system we would expect 
that typical relaxation times are directly related to the 
energy scales that appear in the spectrum, such as the 
bandwidth or the Mott gap. In this case all information 
on the initial energy e is lost for pulses which are short 
enough to resolve the electronic dynamics. The equilib- 
rium interpretation of conventional photoemission data 
in terms of the electronic spectrum of the solid thus be- 
comes meaningless in this limit. 

A full theory for time-resolved photoemission spec- 
troscopy (TRPES), which covers both the nonequilib- 
rium effects of the electronic state and also the conse- 
quences of the frequency uncertainty of the pulse, was 



2 



presented recently by Freericks et al.J^ Their approach 
extends existing theories of conventional photoemission 
spectroscopy to the case where the sample is not in equi- 
librium and measurement pulses have a finite time du- 
ration. The photoemission intensity as a function of the 
probe pulse delay time is related to electronic one-particle 
real-time Green functions of the sample, 19 which fully 
incorporate the noncquilibrium many-body dynamics af- 
ter the pump pulse. This is in contrast to earlier Green 
function approaches^ which treat pump and probe on 
the same (perturbative) level. The relation to real-time 
Green functions allows to make direct contact to recent 
progress in noncquilibrium many-body theory, such as 
the extension of dynamical mean-field theory 2 ^ (DMFT) 
to noncquilibrium . 22 i 23 ' 24 i 25 ' 26 DMFT, which is exact in 
the limit of infinite dimensions^ can provide insights 
into the real-time evolution of strongly correlated sys- 
tems in a nonperturbative way. 

The new one-particle description of TRPES given in 
Ref. [l^ leads to the question whether real-time Green 
functions can be recovered fully from the time-dependent 
photoemission intensity, or whether parts of the elec- 
tronic time evolution are not accessible by TRPES at 
all. Freericks et al. discussed the case when elec- 
tronic equilibrium states are probed by pulses of fi- 
nite time duratiom 12 i 19 This analysis covers experiments 
(e.g., those of Ref. 0) in which changes of external pa- 
rameters such as the electronic temperature determine 
the dynamics of the electronic state, but the probe pulses 
are not short enough to resolve the thermalization of the 
electronic system in response to the pump pulse. For this 
case the electronic state is characterized by its frequency- 
dependent spectrum, and the photoemission intensity is 
given by this spectrum, broadened in accordance with 
the frequency-time uncertainty^ Such a broadening can 
hamper the determination of the electronic spectrum 
from photoemission data; for the experiment of Ref. 0, 
however, it plays a minor role.— 

By contrast, in this paper we investigate TRPES with 
ultrashort pulses that do resolve the thermalization of the 
electrons after the pump pulse. We consider systems with 
a purely Hamiltonian time evolution involving only elec- 
tronic degrees of freedom. The electronic state is then no 
longer characterized only by a frequency-dependent spec- 
trum, but rather by real-time Green functions depending 
on two time variables. We will show from the general the- 
ory of Ref. [l9l that in this case the full time dependence on 
both time variables cannot be recovered from the time- 
dependent photoemission intensity, no matter how the 
pulse length of the probe pulse is chosen. While time- 
resolved photoemission data can be predicted from cal- 
culated nonequilibrium Green functions^ 9 - the converse 
procedure is thus impossible due to the frequency-time 
uncertainty relation. We note that an analogous lim- 
itation is absent in time-resolved optical spectroscopy, 
where a two-time optical conductivity o(t, t') can be mea- 
sured precisely by making the probe pulses sufficiently 
short. In noncquilibrium DMFT, o(t,t') is directly re- 



lated to momentum-averaged real-time Green functions 
under certain conditions^ 

Below we employ the Falicov-Kimball model^ which 
describes localized and mobile electrons on a lattice inter- 
acting via a local Hubbard interaction, to study the rela- 
tion between nonequilibrium Green functions and time- 
resolved photoemission data in detail. We consider an 
idealized setup in which the system is suddenly driven 
out of a metallic or insulating equilibrium state, and 
subsequently relaxes to a new phase due to the Hamil- 
tonian dynamics of the electrons. This model situation 
was recently solved with nonequilibrium DMFTi 2 ^ We 
then study hypothetical time-resolved photoemission ex- 
periments during this relaxation process, and find that 
some aspects of the formation of the metallic state are in- 
deed obscured in the photoemission spectrum due to the 
frequency-time uncertainty. On the other hand, the re- 
laxation of an excited Mott insulator leads to characteris- 
tic collapse-and-revival oscillations^ 9 - which result in os- 
cillating mid-gap weight in the time-resolved photoemis- 
sion spectrum. The above-mentioned uncertainty limita- 
tions notwithstanding, TRPES with ultrashort pulses is 
well-suited to characterize noncquilibrium states of cor- 
related electron systems. However, it will often be nec- 
essary to analyze in detail how the time evolution of the 
Green function translates into the photoemission signal. 

The outline of the paper is as follows. In Sec. [Til we 
briefly outline the microscopic formulation of TRPES 
derived in Ref. [l^ and further discuss the role of the 
frequency-time uncertainty in this theory. We then in- 
troduce the Falicov-Kimball model fSec. IIII[) and discuss 
hypothetical time-resolved photoemission measurements 
on systems that relax to a metallic state fSec.HVp and to 
an insulating state (Sec. [Vj) . The discussion in Sec. IVTI 
concludes the presentation. 



II. TIME-RESOLVED PHOTOEMISSION 
SPECTROSCOPY 

In photoemission experiments with both temporal and 
angular resolution the sample is probed with a finite pulse 
of definite wave vector— q. The detector collects the pho- 
toelectrons which are emitted in a certain direction k ei 
and it is sensitive to their kinetic energy E = fi. 2 fc 2 /2m, 
but not to their arrival time (k e — k e k e is the photo- 
electron momentum). The time- resolved photoemission 
signal is thus proportional to the total number of elec- 
trons per solid angle dfij, and energy interval dE 

T(h - dN(k e ,E;q,t p ) 

I(k e ,E,q,t p )- (1) 

that are emitted in response to a pulse that hits the 
sample at time t p r^ This definition includes only pho- 
toelectrons excited by the probe pulse and omits direct 
photoemission due to the pump pulse. 
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In Ref. an expression for the photoemission sig- 
nal ([TJ) was derived, using only the so-called sudden 
approximation,— which neglects the interaction of photo- 
electrons with the remaining sample. The photoemission 
signal is then only related to matrix elements M(k, q; k e ) 
which couple Bloch states with quasi-momentum k in the 
solid and one-electron scattering states with asymptotic 
momentum k e via absorption of a photon with momen- 
tum q, and to the real-time one-particle Green function 

GlA^')^^[p^ k ,{t')c k {t)]. (2) 

The latter incorporates the full noncquilibrium dynamics 
of the sample: cjj, (i) = U(t,t min yc^U(t,t min ) are an- 
nihilation (creation) operators for electrons in the solid 
with momentum k, whose propagation in time, with 
U(t,t min ) = T T exp[— i dTH(r)/h}, includes all ex- 
ternal fields except for the probe. The initial state at 
some early time t m i n is usually given by the thermal en- 
semble at temperature T, po oc exp[— H(t m i n )/T}. The 
presence of the surface and the dependence of photoemis- 
sion spectra on matrix elements can substantially com- 
plicate the comparison of theoretical and experimental 
data for specific materials. In order to reveal general as- 
pects of TRPES we thus resort to further approximations 
that are commonly made in this context: (i) We assume 
that photoemission measures the bulk properties of the 
sample which are contained in the momentum-diagonal 
Green function G^(t,t') = G^ k (t,t') of the infinite and 
translationally invariant system, and (ii), we take matrix 
elements to be constant but satisfying momentum con- 
servation in the plane, M(q,k;k e ) = MSk n +q h k cir The 
time-resolved photoemission spectrum (JT]) is then given 

I(k e ,E;q,t p ) oc ^ 5 kli +q h k e ^kv(E - cq - $;t p ), (3) 

ka 

T feCT (w;i p ) = -i J dtj dt'S{t)S{t') 

x e^'-^G^it + t^t' + tp), (4) 

where S(t) is the (real) pulse envelope function (centered 
at r = 0) and $ is the work function of the solid. Note 
that Eqn. ([3]) and ([4]) become exact for perfectly layered 
(two-dimensional) structures, when k denotes the two- 
dimensional momentum of the sample. In the following 
we will discuss the momentum- and frequency-dependent 
expression (|4|). For simplicity we will refer to Eq. (J4| as 
the photoemission intensity; observations made for this 
function presumably persist after summation over some 
part of the Brillouin zone [Eq. ©]. 

Eq. simplifies when the system is in equilibrium. In 
this case Green functions depend on the time difference 
only, and the Fourier transform is given by22 

9%M = / dte^G<Jt,0) = 2^L fcCT H/H, (5) 



where A ka (iS) is the equilibrium spectral function^ and 
f{uj) = l/(e w / T +l) is the Fermi function for temperature 
T. The photoemission intensity then reduces to 

Ik*(u) = J ' du'\S{u + u>')\ 2 A ka {J)f{J), (6) 

which is a convolution of the well-known expression 
Zfcer(w) oc Aha (to) f (u>) for the intrinsic photocurrent 
in continuous beam experiment s 17 ' 18 with the Fourier 
transform S(oj) = JdtS(t)e tuJt of the pulse envelope. 
Due to the frequency-time uncertainty of the pulse, the 
frequency-dependent spectrum is thus smoothencd on a 
scale 8ui > 1/St when the pulse has a finite length St, as 
discussed in Refs. [lj^ and[l2l 

Here we study the case of an electronic system that 
is not in equilibrium. The Green function G ka (t,t + s) 
then contains important information both in the absolute 
time t and in the time difference s between addition and 
removal of an electron. However, when the probe pulse 
extends only over a finite length S, the product S(i)S(t') 
in Eq. (0| vanishes for all t — t' > S, and hence G ka (t,t + 
s) enters Eq. (J4| only for s < S. It is therefore impossible 
to deduce G ka (t,t + s) from spectra that were recorded 
with pulses of length S < s. In other words, the time 
resolution (in t) with which t + s) can be measured 

is limited by s. This also becomes clear when attempting 
to invert the convolution of G^ a (t, t') in Eq. (j4}: Starting 
from the Fourier transform 

iha{s\t p ) = J duje luls l k<J (uj;t p ) , (7) 
and using, e.g., Gaussian pulses, 

S(t) = exp(-^), (8) 

we obtain 




While the integral in © apparently measures G ka (t p + 
s/2, t p — s/2) with a time resolution of S, it is practically 
impossible to choose i < s because then the result van- 
ishes compared to any noise added to Xfc CT (s;i p ), due to 
the Gaussian prefactor (whose form is due to © but the 
suppression of the signal for 8 <C s is independent of the 
pulse shape). We conclude that the noncquilibrium two- 
time Green function cannot be fully measured by means 
of TRPES, and one must always carefully analyze how 
the time evolution of the Green function translates into 
the photoemission signal for a given theoretical model. 
This will be illustrated for the Falicov-Kimball model be- 
low. 
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III. THE FALICOV-KIMBALL MODEL IN 
NONEQUILIBRIUM 

In the remaining part of this paper we concentrate on 
one specific model for electronic dynamics in a single 
band, the Falicov-Kimball models This lattice model 
describes itinerant (J,) and immobile (|) electrons which 
interact via the local Coulomb repulsion U. The Hamil- 
tonian is given by 

h = v~ij4i c ji + u ^2 n ii n n - ^ n ^' ( 10 ) 

ij i i<7 

where c^J are annihilation (creation) operators for the 
two species of fermions on lattice site i, and m a = c\ G c iG 
is the corresponding density (a = I, 1). Hopping between 
sites i and j (with amplitude Vij ) is possible only for the 
mobile Q) particles. The Falicov-Kimball model has been 
an important benchmark for the development of DMFT 
in equilibrium, because the effective single-site problem 
for the mobile particles is quadratic and can be solved 
exactly^ 3 - This model currently plays a similar role for 
nonequilibrium DMFT ) 22 i 23 i 24 i 25 i 26 in particular since no 
appropriate real-time impurity solver is yet available for 
the Hubbard model. 

In spite of its apparent simplicity the Falicov-Kimball 
model has a rich equilibrium phase diagram containing 
metallic, insulating, and charge-ordered phases^ In the 
following we consider only the homogeneous phase at 
half-filling for both particle species (nj. = n-f = 1/2), 
which in equilibrium undergoes a metal-insulator tran- 
sition at a critical interaction U = U c on the order of 
the bandwidth . 33 i 34 ' 35 i 36 This phase is studied in an ide- 
alized nonequilibrium situation, by preparing the system 
in thermal equilibrium for times t < 0, and changing the 
interaction parameter U abruptly at t = Oj^ 3 - Of course 
within a more realistic description of the pump pulse the 
system would not be excited into a state that is an equi- 
librium state of any simple Hamiltonian. Nonetheless, 
the interaction quench we study here allows us to dis- 
cuss time-resolved photoemission signals for a situation 
in which the electronic dynamics drives the system be- 
tween different phases. Since the exact DMFT solution 
for the interaction quench is available^ 3 , we can directly 
relate Green functions and photocmisson signals. In par- 
ticular, we will focus on two specific phenomena, namely 
(i) the formation of narrow quasiparticle resonances dur- 
ing the buildup of the metallic state (Sec. IIVj) . and (ii) 
coherent collective oscillations after an excitation of the 
insulating phase (Sec. N} . 

In the following we analyze the time-resolved photo- 
emission signal of the mobile Q) electrons using Gaussian 
pulses ((SJ) and omit the time-independent contribution of 
the immobile (|) electrons. As in Ref. we assume a 
semielliptic density of states p(e) = V '4V 2 — e 2 /(2irV 2 ) 
with half-bandwidth 2V for the single-particle energies 
which are the eigenvalues of the hopping matrix 
Vij. Exact expressions for the real-time Green func- 
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FIG. 1: The momentum-dependent Green function g£ ± (u>) 
[Eq. in the initial (a) and final (b) state for the quench 

from U — 3 to U = 0.5 (n c = nf = 1/2; temperature T = 0). 
Due to the local self-energy in DMFT, G^(t,t') depends on 
k only via e*, in the homogeneous phase, (c) Density nk(t) 
for momentum fe with = — f [thick red curve in (a) and 
(b)]. The horizontal line is at rife (oo). (d) Green function 
GjJ (t + s,t) for the same Differences between the Green 
functions for t = and t — oo are best visible around s = 0; 
their decay is almost identical. 

tions Gu(t,t') of the mobile electrons are given in the 
Appendix (from now on we suppress the index J, of the 
mobile electrons) . We take V = f as the unit of energy, 
so that the full bandwidth is 4 and the critical interac- 
tion is given by U c = 2V = 2f2£ We will also set h = I , 
setting the unit of time as H/V. For example, for V = 
f eV we have H/V = 0.66 fs. 



IV. PUMPING THE INSULATOR INTO A 
METALLIC STATE 

In this section we investigate the formation of a metal- 
lic state in real time, similar to pump-induced insulator- 
to-metal transitions on ultrashort time scales4^^ In the 
Falicov-Kimball model, such a process takes place after a 
quench from an insulating state to the metallic parameter 
regime, which we consider now. We prepare the initial 
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state at U = 3 and temperature T = 0, and change the 
interaction abruptly to U = 0.5 at time t = 0. Subse- 
quently the system relaxes to a new stationary state, in 
which Green functions depend on time difference only^ 
In the following we first discuss the real-time Green func- 
tions for this process and then the corresponding time- 
resolved photoemission signal. 

A. Real-time Green functions 

The difference between the initial and final state is 
evident from the momentum-diagonal Green function 
G"j;(t,t') of the mobile particles [Eq. (|5J)] and its Fourier 
transform 

g< T (Lu) ee t Hrn^ J ds e^ s G< (t + s, t) (f f ) 

in the limit t = ^foo, respectively: While g^_(uj) has a 
broad maximum (Fig. |T^u), a sharp peak in g^ + (uj) indi- 
cates that quasiparticle excitations have a long lifetime 
in the final state (Fig. [TJa). Note that for a quench in 
the Falicov-Kimball model the final state always retains 
memory on the initial configuration^ This memory is 
contained in a nonuniversal occupation function F(uj) 
which replaces the Fermi function f(ui) in Eq. (O, i-e., 
g^ + (ui) = 2iriAk(uj)F(tjj) (see Appendix). 

The development of the metallic state with its sharp 
quasiparticle-like resonances can be observed from the 
full time dependence of the momentum-diagonal Green 
function G^(t,t'). In particular, we characterize the 
transition by means of (i) the total spectral weight, i.e., 
the momentum occupation n k (t) = —iG^(t,t), and (ii) 
the decay of G^{t,t + s) as a function of ,s > 0. From 
the latter one can read off the lifetime of a hole which 
is created at time t in the nonequilibrium state. The 
time evolution of these two quantities is similar for all k; 
it is shown for one representative value of k in Fig. [TJ; 
and[TJi (namely = —1, marked by the thick red line in 
Fig. QJ, and QJ>) : (i) Relaxation of the momentum occu- 
pation rik(t) takes place on a time scale on the order of 
the inverse bandwidth; after this very short time inter- 
val the final value rifc(oo) is almost reached, and a slower 
relaxation follows (Fig. [ft). A similar behavior was ob- 
served previously for the time dependence of the number 
of doubly occupied sites^ (ii) As a function of s, the 
Green function G^(t, t + s) decays slow compared to the 
inverse bandwidth (Fig. QJl). For t = oo, this is in accor- 
dance with the sharp peak in g^ + (uj). However, the fast 
decay is observed even for holes that are created at t = 
+ , when the system is still in the insulating state, indi- 
cating that the lifetime of hole excitations depends only 
weakly on the time of their creation. The reason for this 
behavior is that in the Falicov-Kimball model scattering 
occurs only between mobile and immobile particles and 
is thus determined by the Hamiltonian and the (initial) 
configuration of the immobile particles. If scattering oc- 
curred between two mobile electron species (e.g., as in 
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FIG. 2: Photoemission signal [Eq. Q] (a.u.) for the same 
situation as the Green functions in Fig. [T] (ek = — 1), using 
Gaussian probe envelopes ((SJ) with 8 = 10 (a), S = 2 (b), and 
S = 0.66 (c). The vertical dashed lines are explained in the 
text. t p and S are in units of h/V = 1. 

the Hubbard model), then we would expect the shape of 
the quasiparticle resonances to depend on the quantum 
state of the mobile particles as well, and therefore would 
expect it to change considerably during the relaxation 
process. 

B. Photoemission spectrum 

As discussed in the introduction, it is a central ques- 
tion whether the time-resolved photoemission spectrum 
Ik(ijj]t p ) [Eq. (0}] contains the same information as the 
Green function G^(t,t'). For the present case, the an- 
swer is no: In the last subsection we saw that holes decay 
at a rate T V/h, which depends only weakly on the 
time t when the hole is created, even for short times 
< t < h/V. To establish this behavior from the pho- 
toemission intensity, however, one would have to mea- 
sure Gk(t, t+ 1/r) with time resolution better than h/V, 
which is impossible according to the discussion at the end 
of Sec. HU 

The effect of the frequency-time uncertainty can be 
seen in detail in the redistribution of spectral weight in 
the photoemission signal Tf.(uj;t p ) as a function of the 
probe time t p (for ek = —1, Fig. [5]). When the system 
is probed in a stationary state (t p = ±oo), the intensity 
is given by g^ ± (uj) folded with the spectrum |,S(w)| 2 = 
2ir6 2 exp(-<5 2 w 2 ) of the probe pulse [cf. Eq. ©]. This 
broadening completely washes out the peak for short 
pulses (S = 0.66 in units of h/V, Fig. [2t) . On the other 
hand, pulses much longer than the inverse bandwidth do 
not resolve the fast relaxation which is essentially com- 
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V. OSCILLATIONS OF AN EXCITED MOTT 
INSULATOR 

In this section we investigate the relaxation dynamics 
of a Mott insulator in which a metallic state has been 
created by the pump pulse. In the Falicov-Kimball model 
this can be simulated in the strongly interacting regime 
(U = 10) by preparing the system in a metallic state (U 
— 1, temperature T = 0) at t = 0. Again we first discuss 
the real-time Green functions for this situation and then 
the time-resolved photoemission signal that corresponds 
to them. 



FIG. 3: (a) Relaxation time r(S) [Eq. lfl3]l. 7] = 0.95] for 
the spectral width in the central peak region of the photo- 
emission spectrum at et = — 1 (region between the vertical 
dashed lines in Fig. [2}, plotted against the pulse length 8. 
(b) Ratio Ws(oo)/Ws( — oo). The horizontal dashed line is at 
nfc(oo)/n fe (-oo). 



plete for t > 2, but rather show a time dependence even 
for t > 2 (5 = 10, Fig.[2K). For intermediate pulse length 
neither time nor frequency dependence is resolved (5 = 
2, Fig. [2b). 

For a quantitative analysis we now consider the weight 
in the central region of the peak in the photoemission 
spectrum, 



Ws(t p )= / du;l k (uj;t p ). (12) 



In the present case we use a = —2.39 and b = —2.17 as in- 
dicated by the vertical dashed lines in Fig. [2j they enclose 
twice the full width at half maximum of a Lorentzian fit 
to the peak in Ak{u>). We then define a time t(S) at 
which the relaxation of Ws(t p ) is essentially complete, 

Ws[t(S)}~W s (~^) 
W tf (oo) - W tf (-oo) 1 J 



A. Real-time Green functions 

In the present case the momentum-dependent Green 
function evolves from a well-defined quasiparticle band at 
t = — oo, which is cut off by the Fermi function (Fig. UK), 
to a gapped spectrum at t = oo (Fig. 0b). Note that in 
the final state spectral weight remains in the upper band 
because the system is strongly excited with respect to 
the insulating ground state at U = 10, and there is no 
coupling to an environment to which this excess energy 
could be passed during the relaxation process. 

As in the previous section we consider a representative 
momentum (e^ = 1), chosen such that n k is small in the 
metallic state and increases after the quench (Fig. 2b). 
Again this relaxation takes place on a very short time 
scale (on the order of the inverse bandwidth), but now 
nk passes through a series of damped oscillations with 
period 2ir/U before reaching its final value. These os- 
cillations are characteristic for the dynamics of a Mott 
insulator which is dominated by a Hubbard-type den- 
sity interaction U^ Ji n^riii. In fact, if the Hamiltonian 
were given only by this interaction term, then the time 
evolution operator cxp(— itU n^riii) would itself be 
27r/C/-pcriodic^ and hence oscillations would occur in all 
nonlocal quantities. These so-called collapse- and-revival 
oscillations were first observed and described in experi- 
ments with ultracold atomic gases^ where the Hamilto- 
nian of the system can be designed in a controlled way. 
We will now discuss the fingerprint of these oscillations 
in the time-resolved photoemission spectrum. 



with r/ close to one; we use r\ = 0.95. In Fig. [3^ we see first 
of all that t{8) is proportional to the pulse length for large 
S, which is due to the fact that long pulses merely aver- 
age the spectrum of the final and initial state. Hence t(S) 
does not yield any information about intrinsic relaxation 
times for 6 > 2. On the other hand, for 5 < 2 the ratio 
Ws(oo)/Wg(— oo) approaches the value n k (oo)/nfc(— oo) 
that one would obtain by integrating over the whole spec- 
trum instead of the peak region alone (Fig. [3b) ; this is 
due to the insufficient energy resolution. Therefore only 
the relaxation time of the whole spectral width, i.e., of 
the momentum occupation rik(t), can be determined re- 
liably. 



B. Photoemission spectrum 

In. Fig. [5] the angular-resolved photoemission spec- 
trum Xfc(w;f p ) is plotted for the same fixed momentum 
[efc = 1, using Gaussian pulses ©]. All features of the 
spectrum except for its total weight, which is propor- 
tional to 7ifc(t), are washed out for short pulses (5 = 0.2, 
Fig-Efc), whereas long pulses show the formation of a gap, 
but cannot resolve the oscillating nature of the state (5 
= 0.66, Fig. [5^). For intermediate pulses, however, both 
the 27r/ [/-periodicity and the gap become visible (5 = 
0.33, Fig. Eb)- 
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-2 -1 -2 2 4 6 8 10 
co [V/h] co [V/h] 



1 2 3 4 5 6 

t [h/V] 



FIG. 4: The momentum-dependent Green function g£ ± (uj) 
[Eq. in the initial (a) and final (b) state for the quench 

from U = 1 to U — 10 (n c = n/ = 1/2; temperature T 
u : (c) Density for momentum k with = 1 [thick 

red curve in (a) and (b)]. The vertical dashed lines are at 
multiples of the fundamental oscillation period 2n/U. The 
horizontal line is rafe(oo). 

Interestingly, the coherent oscillations are most pro- 
nounced in the center of the gap (Fig. [5)d). This obser- 
vation can be understood from the atomic limit of the 
Hamiltonian (fTU|) . i.e., for V%j = 0. In the strongly inter- 
acting regime U ^ V the atomic limit gives a good de- 
scription of the transient behavior at short times t < H/V. 
For the interaction term alone, Hjj = U^.n^riii — 
Sio- Mo-T^io- the time evolution of annihilation operators is 
given by e iHut c ja e- iHut = e v -*[c> + (e~ itu -l)c j(T n id -]. 
For t,t' > 0, the Green function then follows as 

G< a (t,t') = ie^-t'Mk* + B ka e u ' u + 

Bl a e- uu + C ka e lU{t '-% (14a) 

with 

A k<J = e ik ^- R i\(l n 3 ,)c) a c la {\ - n lS )) (14b) 

ij 

fl to = Ee* M (»i40o (we) 

ij 

C ka =Y.e lk[R ^ R3 Hn ]S c) (y c l!J n lS ) , (14d) 

ij 

and (■)□ is the expectation value in the (arbitrary) state 




FIG. 5: Quench from U = 1 to U = 10: Photoemission signal 
(a.u.) [Eq. @] for tk = 1, and Gaussian probe envelopes (|8]) 
with 5 = 0.66 (a), S = 0.33 (b), and S = 0.2 (c). Pulse lengths 
8 are in units of h/V = 1. 

at t = immediately after the pump. Inserting this ex- 
pression into Eq. Q we find, for t p 3> 5, 

l k<T (uj;t p ) oc A ka \S{uj + ^ a )\ 2 +C ka \S{uj + n,j - U)\ 2 
+ 2Rc[S(uj)S(oj + fi a - U)B ka e lt " u }. (15) 

The first two terms are centered in the upper and lower 
Hubbard bands at lu = — fi a and u> = U — \i a and do 
not change with time. The third term, which oscillates 
with period 2n/U, has its maximum where S(tu + /i CT ) 
and S(u> + \i a — U) overlap. For Gaussian pulses, this is 
precisely the center of the gap because S(w)S{uj — U) oc 

e -(u-U/2) 2 /2S 2 e -U 2 /5\ 

This discussion shows that on short time scales the ob- 
served time-dependent spectrum in the Falicov-Kimball 
model for large interactions resembled that of the atomic 
limit. Note that the initial state at t — determines 
only the weight of the three components, but not the fre- 
quency of the oscillations. The oscillating midgap weight 
is thus a universal property of the Mott insulator, which 
is largely independent of the excitation process. Via this 
universal feature it may eventually become possible to 
observe collapse-and-revival oscillations in TRPES ex- 
periments on correlated materials. 

VI. CONCLUSION 

In this paper we analyzed hypothetical time-resolved 
photoemission experiments on correlated electron sys- 
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terns that are not in equilibrium, building on the general 
theory of Ref. [H. We showed that the two-time Green 
function, which characterizes the noncquilibrium state of 
the electrons, cannot be fully measured with TRPES, 
no matter how long or short the probe pulses are cho- 
sen. For example, using DMFT for the Falicov-Kimball 
model we found that in the buildup of the metallic state 
the Green function of the transient state cannot be de- 
termined from the photoemission signal. On the other 
hand, if an excited Mott insulator is created by the pump 
pulse, its characteristic collapse-and-revival oscillations 
can nevertheless be inferred because they correspond to 
oscillating weight in the center of the Mott gap in the 
time-dependent photoemission spectrum. 

TRPES is in some sense complementary to time- 
resolved optical spectroscopy, which measures the two- 
time optical conductivity a(t,t'). Under certain condi- 
tions, the latter is obtained in DMFT from a momentum- 
averaged product of two Green functions that also enter 
the expression for the photoemission spectrum^ The de- 
pendence of a(t, t') on both t and t' can be measured 
precisely with sufficiently short probe pulses, unaffected 
by any minimum uncertainty, but unlike in photoemis- 
sion spectroscopy there is no sensitivity toward specific 
momenta k. 

In conclusion, we showed that in spite of the frequency- 
time uncertainty of the probe pulse, TRPES has the po- 
tential to discover fascinating details of the electronic 
thermalization process. Unlike for conventional photo- 
emission on systems in equilibrium, however, the time- 
dependent photoemission signal does not yield the real- 
time Green function directly, so that more detailed com- 
parisons to theoretical predictions are needed. 
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APPENDIX A: CALCULATION OF DMFT 
GREEN FUNCTIONS 

1. Contour Green functions 

In this appendix, which is a direct extension of the 
work presented in Ref. l23l . we give the detailed deriva- 
tion of the real-time Green functions for the interaction 
quench in the Falicov-Kimball model, using DMFT for 
noncquilibriumi^. 

The retarded, advanced, and lesser Green functions are 



defined by 

GL(M') = -m-t')Tr[p {c k<T (t),<iAt')}} ( Ala ) 

GL(M') =i®{t' -t)Tr[ Po {c kl7 (t),c{ a (t')}] (Alb) 

G&r(*,f)= iTr[PocUt')c klT (t)] (Ale) 

respectively, where c ka (t) and po arc defined below 
Eq. @. DMFT for noncquilibrium is based on the 
Keldysh formalism ; 37 ' 38 which yields the contour-ordered 
Green function G ka (t,t') = -i(T c c ka (t)cl a {t')) with 
time arguments on the contour C that runs from t min to 
tmax on the real axis, then from £ max to i m i n , and finally to 
tmin — i(3 on the imaginary time axis {(3: inverse temper- 
ature). Retarded, advanced, and lesser Green functions 
(|A1[) are obtained from the real-time components of the 
contour Green function^ 

GL(M') = [G 1 k 1 a (t,t')-GZ(t,t')] (A2a) 
GL (*,*') - [Gil(t,t') - GH(t,t')] (A2b) 
Gt(t,t') = G k l(t,ti), (A2c) 

where superscripts refer to the two time arguments: 1,2, 
and 3 indicates whether a time argument is on the upper, 
lower, or vertical part of the contour, respectively. These 
relations and also the symmetries 

Gl a (t,t') = GlAt',ty (A3a) 
G< <x (t,t') = -G< a (t',ty. (A3b) 

hold for all contour Green functions considered here. Fur- 
thermore, the contour Green functions obey antipcriodic 
boundary conditions in both contour arguments: 

GH(t, t min ) = -Gll{t, t min -ip), (A4a) 
Gt(t min ,t') = -Gt(t min -0,t'), (A4b) 

for v = 1,2,3. 

For the lattice Hamiltonian (TIT))) , the interacting con- 
tour Green function satisfies a Dyson equation^ 

[{Q k l - S feCT ) *G ka }(t,t') = 6 c (t,t'), (A5a) 

where Sfc (T (t, t 1 ) is the contour self-energy, and G ka (t,t') 
is the noninteracting Green function, whose inverse 

^(MO = 5 C {t,t')[{i$ + AV) - e k(7 (t)] (A5b) 

can be written as differential operator on the contour. 
Here we assumed a homogeneous state, and e ka = 
Y^j V? exp[ik(Rj — Ri)] are the single-particle band en- 
ergies. We also introduced the convolution (f*g)(t,t') 
= J c dif(t,t)g(i,t') of two functions along the con- 
tour, the contour delta function S c (t, t') [defined by 
f c di f(t)5 c (t, t) = f(t)}, and the contour derivative <9 t c .— 
Taking into account the boundary conditions (|A4[) . the 
integro-diffcrcntial equation (|A5[) has a unique solution 
for G krT {t,t'). 
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2. DMFT equations 

In DMFT the self-energy is assumed to be local in 
space, and hence independent of k for a homogeneous 
state. This approximation is exact in the limit of infi- 
nite spatial dimensions^ both for equilibrium and for 
the Keldysh self-energy 32 The local self-energy S (T (t, t') 
and the local Green function G a (t,t') = Gu a (t,t') = 
J2k Gka(t, t') are then calculated from an auxiliary prob- 
lem in which the degrees of freedom at a single lattice site 
i are coupled to some unknown environment, which must 
be determined self-consistcntly. For the Falicov-Kimball 
model, the auxiliary problem is quadrati o 22 i 33 such that 
the equations of motion can be solved explicitly. 

The DMFT equations for the interaction quench in the 
Falicov-Kimball model were derived in Ref. [23|. In the fol- 
lowing we state these equations without derivation, and 
then give details of the solution. In particular we calcu- 
late the momentum-dependent Green function Gk(t,t') 
= Gki(t,t') of the mobile (|) particles, which is needed 
for the photocmission intensity ((¥]). From now on we con- 
sider only properties of the mobile particles, and omit the 
index L 

The local Green function of the mobile particles is 
given by the sum 

G(t, t') = w Q{t 7 1') + Wl R{t, t') , (A6a) 

of the local Green functions Q(t,t') and R(t,t') at sites 
with zero and one immobile (|) particle, respectively, 
weighted with the average density w\ = 1 — Wo of im- 
mobile particles. The functions Q(t,t') and R(t,t') obey 
the equations of motion 

[iflf + n)Q{t,t') - (A * Q)(t,t') = 5%t') , (A6b) 

[iflf + n - U(t)]R(t, f) - (A * R)(t, t 1 ) = 5 c (t, £'), (A6c) 

and boundary conditions (|A4[) . For a quench the inter- 
action is piecewise constant in time, U(t) = Q(t)U+ + 
<d(—t)U-. The effective medium propagator A(t, if) must 
be determined sclf-consistently. For a semielliptic den- 
sity of states of the mobile particles, which we adopt in 
the following, the self-consistency cycle can be condensed 
into closed form ^ 

A(t,if) = V 2 G(t,if), (A7) 

where 2V is the half bandwidth of the density of states. 
Eqn. (|A6[) and (|A7|) form a complete set of equations 
for the local Green function. The local self-energy E 
of the mobile particles is then obtained from the Dyson 
equation of the local problem, [(i<9f + n)G(t, t') — [(A + 
E) * G](t,t') = S c {t,t'). Together with Eq. ifAl)]) . this is 
easily transformed into 

wiU(t)R(t,t') = (E*G)(t,t'). (A8) 

Finally the local self-energy £(£, t') is inserted into the 
lattice Dyson equation (|A5[) . 

(iflf + e k )G k (t, t') - [E * G k ] (t, t') = S c (t, 0, 

(A9) 



which yields the fc-dependent Green functions. Note that 
t') depends on momentum only via the single parti- 
cle energy because we assumed a homogeneous state. 

3. Langreth rules 

To solve the contour equations (|A6b[) , (|A6c[) , (|A8p , and 
(|A9[) we first rewrite them in terms of their retarded and 
lesser components, using the identities (|A2[) . In effect, 
this means that contour derivatives C(t,t') = d%A(t,if) 
are replaced by real-time derivatives,— 

C I (t,t') = d t A I (t,t') (AlOa) 
C < {t,t') = d t A < {t,t'), (AlOb) 

and convolutions C(t,t') = [A * B](t,t') of two contour 
Green functions A and B are expressed in terms of their 
retarded, advanced, and lesser components according to 
the Langreth rules 3 ^ 

C x (t,t') = J dfA r (t,£)B r (t,i'); (Alia) 

C < {t 1 i') = f diA < (t,t)B a (i,t') 

+ f diA T (t,t)B < (i,t'). (Allb) 

•J — OO 

The integral boundaries account for the fact that re- 
tarded (advanced) Green functions A 1 ^ (t,t r ) vanish 
when t < if (t > t'). Furthermore, we shifted £ m j n — > 
—oo in the second equation, such that the convolution 
extends over the whole axis but contributions from the 
vertical part at i m i n — it can be dropped. This step is 
discussed in further detail below. 

The contour delta function on the right-hand-side of 
(| A6bp . (|A6c|) , and (|A9[) vanishes when the lesser com- 
ponent is taken, and it is replaced by the usual delta 
function S(t — t') for the retarded components. However, 
because any retarded function A r (t, t') vanishes for t < 

retarded equations of motion are only considered for i 
> if ', and the initial value at t — t' is determined by the 
weight of the delta function and the derivative operator. 
In particular, we obtain 

Gl(t, i) = R r (t, t) = Q T (t, t) = -i (A12) 

from Eqn. (|A"6b"|l . (|A"ocJ) . and (K9\i . These conditions 
follow also directly from the anticommutation relation of 
creation and annihilation operators. 

4. Stationary states 

For the interaction quench we treat the equations of 
motion separately in the four regions where both t and 
if do not change sign; we introduce additional subscripts 
+ and — which indicate whether the time arguments are 
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greater or less than zero, respectively. Inserting (|A10aP 
and (|Alla[) into Eq. (|A6[) yields a closed set of equations 
for R T (t,t') and Q r (M'), 

A T (t,t') = V 2 [ Wl R r (t,t') + w Q r {t,t')\ (Al3a) 
[idt +fj]Q T (t,f) = f dtA z (t,t)Q T (i,t') (A13b) 



[id t + n-U^R 1 ^^') = J dtk T {t,^)R T {t,t'), (A13c) 

which must be solved for t > t! , using the initial condi- 
tion (|A12[) . The self-consistency equation (|A7|) was used 
in Eq. (|A13ap . Note that in Eq. (|A13|) . Green functions 
with both time arguments greater or less that zero, i.e., 

the (-) — h) and ( ) components, do not mix with other 

components. Because U(t) is constant for t > and £ 
< 0, respectively, the solutions of (|A13[) are thus trans- 
lationally invariant in time when both t and t' have the 
same sign, and we make the ansatz 



A* ±± (t,t')=a* ± {t~t') 



(Al4a) 



pOQ 

a T ±(z)= ds e lzs a\(s) (A14b) 
Jo 

for all contour functions A — G, R, Q, A, Gk, and £ 
(with a = g, r, q, A, gk, and a, respectively). Using this 
ansatz in Eq. (|Al3j) we obtain a set of cubic equations, 

g±(z) = w q±(z) + wiq±(z) , (A15a) 
^(zHIz + m-U 2 .^)]- 1 (A15b) 
f%{z) = [z + f i- V 2 g* ± {z) - C/i]- 1 . (A15c) 

that can be solved analytically. These cubic equations 
are well-known from the DMFT solution of the Falicov- 
Kimball model in equilibrium.— This is of course ex- 
pected when both t and t' < 0, because before the quench 
the system indeed is in an equilibrium state. In a similar 

way, the retarded (-1 — h) and ( ) components of £ and 

G k are obtained from Eq. (|A"gj) and (|Al)]) . 



(A16) 
(A17) 



<5±(z) = wiU±r±(z)/g±(z) 
9k±( z ) = [z + fi-e k - cr^(z)] -1 . 



Furthermore, advanced Green functions are directly re- 
lated to the retarded ones by symmetry (|A3j) , so that we 
have 



Al ± (t,t') = al(t-t') 



a±(z) 



dse izs a a ± {s) = d T ± (z*y 



(Al8a) 
(Al8b) 



The lesser Green functions are translationally invariant 

in time only for both t and t' < [( ) component], 

when the system is still in equilibrium. One then has2£ 

A<(t,t') = [ ^ e ^*'^a<H, (Al9a) 
J 2n 

a<{cu) = f(w)[al (w) - a r _{uj)}, (A19b) 



where f(w) = l/(e"/ T + 1) is the Fermi function, and 
— a T _(u>) = — 2iImaL(w) is proportional to the 
spectrum of the equilibrium Green function. Mathemat- 
ically this follows from the solutions of the equations of 
motion on the full contour, including the vertical part, 
and taking into account the antiperiodic boundary con- 
ditions. For the quench we use (|Al9j) as initial condition 
for the lesser components; only then can we then let i m ; n 
— > — oo, and disregard the vertical part in the Langrcth 
rule (pnibjl . 

On the other hand, we show below that in the limit 
where both t and t 1 tend to oo (but their difference is 
finite), the lesser (++) components take a form very sim- 
ilar to (|A19|1 . 



lim Af , (t + s,t) = 

t—>oo 



du> 
2^ 



(A20a) 



a<{uj) = F(w)[a%(u) - a T + (uj)]. (A20b) 

The function F (u>) is common for all a = g, r, q, o~, and 
gk ■ One can in fact directly see from the equations of mo- 
tion (|A6bj) . (|A6cj) . l[A"8"j) . and (|A9|) for the lesser compo- 
nent that if the stationary limit (|A20a[) exists, then Green 
functions a+(w) must have this common factor F(ui). To 
find this factor, however, the equations of motion must be 
solved, because it contains the entire information about 
the initial state. 



5. Double Fourier transforms 

We now consider the cases with one or two positive 
time arguments, i.e., after the quench. We introduce 
double Fourier transforms 



A\_(z,n)= / dt 



die 1 



dt' e 1 "*' A r + _(t,t') (A21a) 

O 

dt' e"'*' 'At+tf ,t), (A21b) 



JO J -oo 

for retarded and advanced components, 

-^e- wt 'A<„(z,Lu)= dte lzt A<_(t,t') (A22a) 
2 71 " Jo 

7 />00 

■^-e.- lut 'A< + {uj,z) =J dte tzt A< + {t',t) (A22b) 

for the lesser components with mixed time arguments 
(which holds for t' < 0), and 

poo poo 

A< + (z,r))= dte lzt \ dt' e tTlt 'A< + (t,t') (A23) 
Jo Jo 

for the lesser Green function with both time arguments 
after the quench. In this subsection we derive explicit ex- 
pressions for A r ^ (z,u)), A+_(z,lo), and A + + (z , u); the 

remaining are then obtained by symmetry (|A3[) . 



£L + {r,,z) = A\_(-if-z*y 
A< + (iv,z) = -A<_(-z*,-wY 



(A24a) 
(A24b) 
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Using Langrcth rules (|A11|) once again yields for convo- 
lutions C = A* B, 



J 



(z, -u) = A 1 , (z, -oj)b T _ (w) + af + (z)B T , (z,-u) 



C<_(z, -cj) = A<_(z, -uj)b a _(u) + a T + (z)B<_(z, -u) + A r + _(z, -u>)b<(u) 
C<4z,r,)=A<4z, V )b%(- V )+V + (z)B< + (z,r,)+ [ ^ [A<_(z, -u)B a _ + (u;, 77) + A\_{z, -u>)B< + (u>, r?)] 



(A25a) 
(A25b) 



2tt 



Furthermore, the derivative C = d%A{t, if) translates into 



C I + _(z,r 1 ) = zA I (z, V )-iaf_(- V ) 



C<_(z,u) = zA<_(z,-u>) 



(A26a) 
(A26b) 



zA< + (z,ri) 



iafUlfj) 
^A< + ( U ,r,), (A26c) 



where one must use the continuity of the components 
at the boundary t = and if = 0. e.g., A r + _(0, t') = 

AL_(0,f)- 

Using (|A25|) and ()A26|) one can rewrite Eq. (|A6j> for 
the various components and solve them using the self- 
consistency equation (|A7[) . For instance, we obtain the 
(H — ) component of the retarded Green functions as 

R r + _(z, V ) = [A l + _(z,7 1 )+ l y + (zy_(-r 1 ), (A27) 
Q r + _(z,r,) = [k\_{z^) + i]q\{z)q*_{-r ] ), (A28) 



where Eq. (|A15[) was used once. Together with the self- 
consistency (|A7j) . this is a simple linear equation for 
A 1 ^ (z,rj). In a similar way all components are deter- 
mined successively: Starting from the retarded (H — ) and 
advanced ( — h) components, the results enter the lesser 
(H — ) and lesser ( — h) components [cf. Eqn. (|A25b[) and 
(IA26bp j. which in turn enter the equations for the lesser 
(++) component [cf. Eqn. (|A25c|) and (|A26c|) ]. The 
procedure is repeated for Eq. (|A6|) . Eq. (|A8|) . and fi- 
nally for the lattice Dyson equation (|A9[) . which yields 
the momentum-dependent Green function Gk(t,t')- 

For completeness we state the final result for G, R, and 
Gfc. For this we introduce the abbreviations 

Kj = [1 - V 2 { Wl f x a rl + w o & )}-\ (A29a) 
< = r x a q x J~9 x a - (A29b) 
Kl = 1 + w lWo U + U. K x aK y M x l (A29c) 

with superscripts x,y <E {r, a}, and subscripts a, /? € 
{+, — }, and we use the convention that the variables of a 
function a x a is (i) z when x — r and a = + , (ii) —77 when 
x = a and a = +, and (iii) u> when a = — . In terms of 



(A25c) 



I 

these expressions the final result is 

G r + _(z,-c) = zU- 2 (A/; r _-l), 



R t + _(z,-uj) = ir T + r T _M TT 



+- > 



G k+ _{z, -u>) = ig\ + gl_K%_ 



and 



(A30a) 
(A30b) 
(A30c) 

G<_(z, -u) = if(u){M™_ - A/; r _)/U 2 , (A31a) 

R<_(z, -u) = if{ujy + {ftM™_ - r x _M"_), (A31b) 

G< + _(z, -oj) = if(u,)gi + (g a k _K™_ - g T k _K^_). 

(A31c) 

The lesser (H — h) component can be written in the form 

A< + (z,r ] ) = F(z,7 1 ) a+ ~ a+ +F A (z,r,) (A32a) 
rj + z 

for A = G, R, and G k , where 



M aa_ 



Ml 



F [Z,V)= / , , T/2 /-a =i 

2tt z + r] + V z \gl. — g. 



AH 



and 



F G (z, ?? ) = F(z,7 7 )/U 2 , 

x (f a j^jf; a _ - f r _A'/; r _Af^ r _), 



(A32b) 
(A32c) 
(A32d) 



o^i^+k+k+^^z+^aT- - K r _A/; r _A.T_) 

+ ~gl- AT- - ft- *T-] } • (A32e) 



6. Back transformation 

To obtain the physical real-time Green functions, we 
have to invert the double Fourier transformations (|A21|) . 
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(|A22|) . and (|A23|) . using the final expressions (|A29|) 
through (j A32|) . Here we give an explicit formula for the 
partially Fourier-transformed lesser component 



A<(uj,t) = / dse ius A < (t + s,t), 



(A33) 



(A = G, R, and Gk)- The singularity at 77 + z = in 
Eq. (|A32a[) determines (|A33[) in the limit t — > oo. 



a+(o;) = lim A < (w,t) 

' t— >oo 



Im#(w,-w)[a+(w) - a+(w)], 



(A34) 



which is of the form discussed above [cf . Eq. (|A20j) ] , with 
F(uj) = ImF(w, —uj). The full result is given by 



A<( w ,t) = £e(±t) 



a±(w) + e 



2vr 



(A35a) 



where 



A_(o;, r?) = i<_(u>, r?) + 5 ~ ( . ^ f H , (A35b) 



A + (uj, 77) = 2ilm 



F(uj, V )[a%(- V )-^ + (uj)] 

UJ + T) 



(A35c) 

27TI UJ — ill — UJ 



and the components for A — G, R, and Gk were given in 
the previous subsection. Note that the first term on the 
right hand side of Eq. (|A35c[) is regular at 77 = —uj, be- 
cause both F(tu, 77) and 55. (—77) — a^_(a;)] are then purely 
imaginary. 
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